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We propose an implementation of a quantum router for microwave photons in a superconducting 
qubit architecture consisting of a transmon qubit, SQUIDs and a nonlinear capacitor. We model 
and analyze the dynamics of operation of the quantum switch using quantum Langevin equations 
in a scattering approach and compute the photon reflection and transmission probabilities. For 
parameters corresponding to up-to-date experimental devices we predict successful operation of the 
router with probabilities above 94%. 
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INTRODUCTION 

Given recent advances in experimental quantum com¬ 
putation, where few-qubit systems have been realized [1- 
6] , a considerable amount of research is currently devoted 
to investigating larger-scale systems such as networks for 
quantum communication [7-9]. In such systems, where 
information must be coherently transported over long 
distances, photons are suitable candidates as quantum in¬ 
formation carriers because of their long coherence times. 
Solid-state devices, on the other hand, seem preferable 
for storage of quantum information [10, 11]. Essential 
components in any quantum communication toolbox are, 
thus, devices capable of directing photons through differ¬ 
ent channels [8, 9, 12-14]. These devices include single¬ 
photon transistors [15-17], switches or routers [18-22], 
single-photon beam splitters [23, 24], etc. 

So far, several approaches have been proposed [25] for 
building a quantum switch, such as an optical implemen¬ 
tation using polarized photons and trapped atoms or a 
phase gate implementation [25]. Also, proposals for a 
beam splitter based on Superconducting Quantum Inter¬ 
ference Devices (SQUIDs) to route photons [23, 24, 26- 
28] or using toroidal resonators [29] have been put for¬ 
ward. 

Here we propose a solid-state implementation of a 
quantum router using superconducting qubits. Given an 
n-photon input, consisting of a train of photon pulses 
where each of them can be in two states (or in a coherent 
superposition of both), the quantum switch absorbs the 
first photon and forwards the next n — 1 photons to a 
path determined by the state of the first absorbed pho¬ 
ton. Our proposed device can be integrated in a larger 
network, where the output of one router is the input of 
the next. As a result, since the routers are controlled by 
the input signals there is no need to control the network 
externally. Also it requires the same number of photons 
to send a signal through two different paths with the 
same number of nodes, while other proposals [9, 18] may 
require less photons for some preferred paths, leading to 
a network where some routes (e.g., a route that at each 


bifurcation node takes the rightmost output path) have 
higher efficiency than others. We analyze the dynamic 
operation of the router using quantum Langevin equa¬ 
tions combined with input/output scattering formalism, 
taking into account decoherence due to relaxation and 
dephasing. We predict successful operation of the router 
with probabilities above 94% under realistic experimen¬ 
tal conditions. 

The paper is organized as follows. First we present our 
proposal for the superconducting circuit that operates as 
a quantum router and derive the effective Hamiltonian of 
the system. We then analyze the dynamics of operation 
of the quantum switch using a scattering approach and 
calculate the probabilities of reflection and transmission 
of an incoming photon. In the last section conclusions 
and a discussion of the possible applications of this 
quantum device are presented. 


MODEL 

A schematic of the device we propose, based on circuit 
quantum electrodynamics (cQED), is depicted in FIG. 1. 
The device is composed of three transmission lines capac- 
itively coupled to four SQUIDs and a transmon, acting 
as artificial atoms that absorb and reflect or transmit the 
photons forward. The device operates as follows: the first 
photon of a register that arrives at the switch is absorbed 
by the transmon which, after being excited, modifies the 
energy spectrum of the SQUIDs in such a way that the 
next photons can only be absorbed by the two SQUIDs 
labeled with 2a or 2b (3a or 3b) if the transmon is in its 
first (higher) excited state. 

The excited SQUID then decays while emitting a pho¬ 
ton into its corresponding outgoing transmission line. 
The transmon also decays emitting a photon into the in¬ 
coming transmission line. The capacitances C 2 sa, G 2 sb, 
Cssa and Cssi, are chosen such that they are smaller than 
any other capacitance in the system. In this way, the 
transmission of the control photon, which is the element 
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He describes the exchange interaction between the 
transmon and SQUIDs with the transmission lines, which 
are modeled as a bath of harmonic oscillators [31, 32] with 
ladder operators 6 i, 62 and 63 for the first, second and 
third transmission lines, respectively: 
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FIG. 1. Circuit QED device proposed to operate as a quantum 
router. This device is composed of four SQUIDs, each with 
a capacitor {C 2 sa, C' 2 s 6 , Cssa, Csab) and a pair of Joseph- 
son junctions (with energies Ej 2 a, etc.). It also contains a 
transmon qubit, consisting of a capacitor (Ct) and a pair of 
Josephson junctions iEjt). These five elements are capaci- 
tively coupled to an incoming transmission line (Ui) and to 
two outgoing transmission lines (V 2 and V 3 ). The fluxes ipq 
in the nodes of the circuit, are quantized variables (see [30] 
for the quantization of these variables) and describe the ab¬ 
sorption and emission of incoming photons by the transmon 
and SQUIDs. 

of the register that controls the operation of the router, 
is prevented [30] . Also the coupling strength between the 
SQUIDs and their farthest transmission line is strongly 
reduced if this condition on the capacitances is satisfied. 

The dynamics of the router in FIG. 1 is described by 
an effective Hamiltonian [30] 

Heff=Hays+He + HT. ( 1 ) 

Here Hsys describes the energy levels of the transmon 
and SQUIDs and also the interaction between them: 
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with k G {2a, 2b, 3a, 3b}. Here a'*' and a denote the cre¬ 
ation and annihilation operators for each of the energy 
levels of the system. These operators create an excitation 
with energy luti in the transmon or uik in the SQUIDs. 
The last term describes the density-density interaction 
between the i-th level of the transmon and the excita¬ 
tions in the fc-th SQUID, with interaction strength Jik- 


The interaction strength of this coupling is given by tti 
and tt 3 , which are the lifetimes of the excited levels of the 
transmon, and also and Tb, which are the lifetimes of 
the SQUIDs. Note that the transmon is only coupled to 
the incoming transmission line and that the second level 
of the transmon is not coupled (directly) to the trans¬ 
mission lines [33]. The coupling between the third level 
of the transmon and the transmission line —which is not 
present in the Jaynes-Cummings Hamiltonian, the ex¬ 
pression that usually describes cQED systems similar to 
FIG. I [2, 4, 34]— is achieved by introducing a nonlinear 
capacitor in the transmon [35] (see also Sup. Mat.). This 
nonlinear capacitor can be realized, e.g., by placing car¬ 
bon nanotubes between the plates of the capacitor [36]. 
These give rise to an energy spectrum that has the form 
EiV) = -b aV‘^) for small voltages U[37, 38]. 

The energy spectrum of the device, described by the 
Hamiltonians Hgys and He [Eq. (2) and (3)], is shown 
in FIG. 2. In this figure the three levels of the trans¬ 
mon are represented with solid lines. The energy levels 
of the SQUIDs, in dashed lines (coupled to the second 
transmission line) and dotted lines (coupled to the third 
transmission line), do not coincide since all the capaci¬ 
tances and Josephson energies are different. After send¬ 
ing a control photon —with energy or ujt 3 — into 
the router (see Fig. 3), due to the coupling between the 
transmon and SQUID energy levels described in Eq. (2), 
the energy needed to create an excitation in SQUID 2a 
when the first level in the transmon is occupied can be 
made equal to the energy needed to create an excita¬ 
tion in SQUID 3a when the third level in the transmon 
is occupied (likewise for 2b and 3h). This can be done 
by tuning the energies of the Josephson junctions of the 
system. Thus, the transmon —the control element of the 
router— forwards a photon with energy uja (or ujb) to the 
second or third transmission line, depending solely on the 
state of the control photon. 

The transmission lines are described by a continuum 
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FIG. 2. Schematic representation of the energy spectrum of 
the Hamiltonian in Eq. (2). This spectrum only contains the 
‘single-photon levels’, i.e., the levels that are accessible only 
by absorbing one single photon (plus |T2), for completeness). 
Energy levels represented with solid lines are coupled only to 
the incoming transmission line. Dashed lines describe levels 
coupled to the second outgoing transmission line and also the 
incoming one. Dotted lines describe levels coupled to the 
third outgoing transmission line and also the incoming one. 
The separation of the energy levels is not to scale. 


We start by deriving the Langevin equations from the 
effective Hamiltonian i?e// [Eq- (1)] including relaxation 
and dephasing [30, 39]. We then obtain the equations 
of motion within the input/output formalism of quan¬ 
tum optics [16, 32, 40]. From these equations we obtain 
the scattering amplitude of transmission and reflection 
of photons, assuming a Lorentzian pulse shape for the 
incoming photons [16]. 

Due to the multiple energy levels of the transmon (de¬ 
picted in FIG. 2), the equation of motion for this element 
of the system —in terms of the operators oi and a\ intro¬ 
duced after quantizing the fluxes [30] — is a complicated 
equation to work with. In order to simplify the calcu¬ 
lations we introduce three projection operators an, an 
and an and their hermitian conjugate —whose defini¬ 
tion can be found in the Supp. Mat.— to replace oi and 
a\. Moreover, we have limited our Hilbert space by con¬ 
sidering only three excited energy levels in the transmon 
and only one in each of the SQUIDs. These operators, 
with the properties 
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FIG. 3. (Color online) Energy spectrum of the Hamiltonian 
in Eq. (2) after the transmon has been excited, (a) After 
the first energy level of the transmon has been excited, the 
next photon, which has energy either uJa or utb (or a super¬ 
position of both) can only excite the energy levels coupled to 
the second outgoing transmission line, (b) If the third energy 
level of the transmon has been excited, the same photon can 
only excite the energy levels coupled to the third outgoing 
transmission line. 


of oscillating modes as [31, 32]: 

Ht= J dpp (b\{p)bi{p) + bl{p)b2{p) + bl{p)b3{p)^ . 

(4) 

ANALYSIS 

In order to analyze the dynamics of the quantum router 
we have studied the scattering of one and two photons by 
our proposed device. Here we calculate the probabilities 
of reflection and transmission through each of the output 
channels. 


and also 

an \Ti) = 1G5) 

an |G^) = 0 

an \Tj) =0 j, (6) 

with IGS”) being the ground state of the transmon and 
\Ti) its excited states, are projection operators that do 
not satisfy the usual commutation relations, so special 
care has to be taken when working with them. On the 
other hand, one of the advantages of using these three 
operators is that, although we will have three equations 
for the transmon instead of one, they are easier to solve. 
The use of these three operators instead of the original 
one does not change the model as long as the number of 
excited states of the transmon is limited to three. An¬ 
other advantage of using this notation is that the final 
quantized Hamiltonian takes a simpler form [the results 
being Eqns. (2) and (3)] and is easier to interpret: the 
operator describing the occupation of the z-th level of the 
transmon is simply atp^an- 

Before deriving the Langevin equations we must find 
an expression that describes the incoming and outgoing 
photons in the input/output formalism. Following [40] 
we introduce the operators b\^^ and ^'^‘4 their her¬ 

mitian conjugates, which create an incoming and an out¬ 
going photon in the second transmission line (the same 
procedure must be applied to the other two transmission 
lines). These operators are defined as (see e.g. [16, 40]): 
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where uj is the frequency of the photon and & 2 (w, ti) is an 
initial (or final) value of the Heisenberg operator defined 
at a time ti —>■ —oo for the input operator (or ti —>■ +oo 
for the output operator). The Langevin equations of the 
‘static’ variables of the system in terms of bun, b 2 in and 
bsin are shown in the Supp. Mat. Once the Langevin 
equations are found, relaxation and dephasing ratios are 
then introduced into this set of equations by following 
Refs. [39, 40]. 

Within the input/output formalism, the scattering am¬ 
plitude for a photon with frequency k in transmission line 
1 to be reflected or transmitted in a transmission line i 
with frequency p is given by 

S{p, k) = (0| bi^out{jp)b\,^{k) |0). ( 8 ) 



ojti 

WtS 

W2a 

bJa 

Refl. even modes 

0.997 

0.947 

1.9 • 10“^ 

1 

Refl. right-moving 

0.998 

0.973 

0.238 

1.6 • 10"® 

Transmission 2 

0 

0 

0.952 

9.8- 10"® 

Transmission 3 

0 

0 

1.6 • 10“® 
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TABLE I. Probabilities of reflection and transmission of an 
incoming photon. The probabilities in this table have been 
computed for different frequencies of the incoming photons 
and considering both even modes and right-moving modes. 
Similar values can be obtained for left-moving modes. The 
probability of reflection of a right-moving photon, as com¬ 
puted here, is equivalent to the probability of absorption plus 
reflection of the incoming photon. 


In the case that multiple photons are sent in, the scat¬ 
tering amplitude reads 


S{pi,...,Pn-,ki,...,kn) = 


— ( 0 | hi^outipi) ■ ■ ■ bi!^out{pn)b\j^^(ki)... b\^^{kn) jo). 

(9) 


These scattering amplitudes are valid for incoming pho¬ 
tons with a definite frequency ki,... ,ki each or a super¬ 
position of these. Assuming that the incoming photons 
have a frequency distribution given by a Lorentzian cen¬ 
tered at ki = LOi and with width l/r^, the scattering 
amplitude of a photon with frequency wi is given by 


/3(p) = 


dki 1 

i{uJi — k) 


-YS{pi,ki). 

Tl 


( 10 ) 


The probability of this process to happen is given by: 


conditions uj 2 a - Ji 2 a = W 3 a - J 33 a = and, similarly, 
^ 2 b — Ji 2 b = bU 3 b — J 33 b = bJa- This cnsures that the router 
is controlled by the state of the first incoming photon. 

The results of our calculations are shown in Tables I 
and II. The first table contains the probabilities of trans¬ 
mission of an incoming photon with different frequencies 
wti, wt 3 , u! 2 a and Wa through the second and third trans¬ 
mission lines together with the probabilities of reflection 
in case this photon is described as an even mode or as a 
right-moving mode (see Sup. Mat. and [32] for a descrip¬ 
tion of these modes). In the case of a photon with energy 
uiTi or lots, since both energy levels are coupled only to 
the incoming transmission lines [see Eq. (3)], we find a 
large probability of reflection, above 94%. It is interest¬ 
ing to notice the different behavior of the router with 
respect to the form of the modes. In case the photons 
are described by even modes, the probability amplitude 
of reflection scales as 


P = J dp\f3{p)f . (11) 

Let us consider again the single photon processes. For 
a good performance of the quantum router, the first pho¬ 
ton to arrive must be absorbed by the transmon and be 
reflected back after all the other photons have been trans¬ 
mitted. It is, thus, necessary to find a procedure to find 
out whether the photon has been reflected after being 
absorbed or without being absorbed. For this purpose 
we describe the incoming photons as even modes given 
by the superposition of left and right-moving modes as 
biin{k) = {riin{k) +hin{k))/y/2 [32]. With this proce¬ 
dure one can see, after computing the scattering ampli¬ 
tudes and the probabilities, that a photon described as 
an even mode is absorbed (and reflected) with a large 
probability (see Tables I and II below). Moreover, the 
probability of transmission of subsequent photons is en¬ 
hanced if these are even modes. 

The probabilities of reflection and transmission of pho¬ 
tons are obtained by numerically integrating the squared 
scattering amplitude over all momenta and imposing the 


1 - f 7TTi 
1 + Itt'ti 


( 12 ) 


for a photon with frequency ujti- In this expression, 7 
is the rate of decoherence and tti is the lifetime of the 
excited energy level in the transmon. In case of right- 
moving modes, since a right-moving mode is reflected 
as a left-moving mode if it is absorbed and as a right- 
moving mode if it is reflected but not absorbed [16], the 
corresponding expression in Table I represents the prob¬ 
ability that a photon is absorbed and reflected, with an 
amplitude that is given by 


bright 


1 + jITti 


(13) 


A photon will only be transmitted if it can excite one 
of the SQUIDS. Table I also shows that a photon with 
energy W 2 a, which coincides with the energy of the ex¬ 
cited state of one of the SQUIDs (see FIG. 2), is totally 
transmitted if it is described as an even mode and only 
partially transmitted if it is a right-moving mode. Since 
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TABLE II. Probabilities of reflection and transmission of a 
photon while one of the levels —Tl or TS— of the transmon 
is excited. Only even modes have been considered (see the 
text). 

we are only interested in a scenario where the photons ab¬ 
sorbed by the SQUIDs are always transmitted, we only 
consider even modes from now on. Notice that the prob¬ 
ability of transmission of a photon indeed goes to one in 
case it is transmitted to the nearest transmission line of 
the excited SQUID (with even modes) and vanishes in 
case it is transmitted to the other outgoing transmission 
line. 

More interesting is the case where two photons (a con¬ 
trol and a target photons) are sent into the router. In 
this case, while the transmon is in the (slow) process of 
absorption and emission of the control photon, the sec¬ 
ond (target) photon is absorbed and emitted in a faster 
process by a SQUID before the transmon returns to its 
ground state. The data corresponding to these scatter¬ 
ing processes is shown in Table II. In this case, a photon 
described by an even mode with energy either uJa or uJt is 
transmitted to the second transmission line with a prob¬ 
ability above 95% if the first level of the transmon is 
excited. The probability amplitude of transmission of a 
photon with energy to the second outgoing transmis¬ 
sion line is then given by 


If the third instead of the first level is excited, the 
photon is forwarded to the third outgoing transmission 
line with the same probabilities as in the case a photon is 
forwarded to the second transmission line when the first 
level of the transmon is excited. For any other frequency 
of the incoming photons these are not transmitted. 


CONCLUSIONS AND OUTLOOK 

We have proposed and analyzed a transmon-based 
quantum router containing a nonlinear capacitor that op¬ 
erates with photons in the microwave regime. The non¬ 
linear capacitor is the element responsible for the photon- 
transmon interaction. We predict successful operation of 
the quantum router with probabilities above 94% for cur¬ 
rent experimental parameters. That is, the probability 


that a photon is transmitted successfully to the target 
channel, assuming that the transmon has been excited 
with a high probability (as the calculations suggest) and 
that the photon is transmitted before the transmon re¬ 
turns to its ground state, is above 94%. This router can 
be built with a set of capacitors and Josephson junctions 
(see Ref. [30] for a proposal of such parameters) in which 
two of the SQUIDs (coupled to the same outgoing trans¬ 
mission line) can be excited by a photon with frequency 
LOa or uJb if one (or the other) level of the transmon is 
excited. 

The quantum router can be used as a single-photon 
transistor, with the distinctive characteristics that both 
the control and target photons can come from the same 
transmission line (i.e., from the same source) and that 
the target photon can be transmitted into two different 
transmission lines. The control photon thus not only 
controls whether the target photon is transmitted or not 
but also into which transmission line it is routed. Since 
the router is operated quantum-mechanically, this allows 
for transmission of photons in a superposition of paths if 
the control photon is in a superposition of states. Because 
of the latter, the proposed router can also be used to 
construct a quantum random access memory [41], where 
a tree-like network with quantum switches at the nodes 
leads an address register of qubits from the root node to a 
superposition of memory cells. Finally, the router could 
be used as a basic element in quantum communication 
networks. A list of requirements to fulfill this purpose is 
discussed in the Supplemental Material [30]. 

This work is part of the research programme of the 
Foundation for Fundamental Research on Matter (FOM), 
which is part of the Netherlands Organisation for Scien¬ 
tific Research (NWO). 
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Proposal for a transmon-based quantum router: Supplemental Material 


These notes contain some comments and additional information needed to understand the derivation of the results 
presented in the main text. We discuss in detail the need for a nonlinear capacitor in the proposed router design 
and possible candidates together with their properties. We also show how our model (Hamiltonian) has been derived 
and quantized and how scattering amplitudes and probabilities are computed. We present the set of numbers we 
have used and discuss obstacles and possible remedies towards scaling up of our proposed router. Finally, we discuss 
the requirements that we believe a quantum router needs to satisfy in order to construct a network for quantum 
communications. 


NONLINEAR CAPACITORS 


In circuit QED, a device similar to the one proposed in the main text (see also FIG. S4) is usually described 
by the Jaynes-Cummings (or Tavis-Cummings) Hamiltonian [1-3]. This Hamiltonian only allows the excitation of 
one level of the transmon by absorbing one photon. If the third level of the transmon is to be reached, then three 
photons are needed. Here we propose a mechanism with which the third level of the transmon can be excited in 
one step by absorbing one photon. Unfortunately this transmon-cavity interaction cannot be achieved by using 
Josephson junctions or any other inductive element (when including an inductive element in the Lagrangian and 
performing a Legendre transformation, the resulting Hamiltonian does not contain a term coupling the transmon 
with the transmission lines, as is required). The nonlinearity provided by a Josephson junction is insufficient to 
couple the third level of the transmon to the transmission lines in a way that allows for a photon to be absorbed 
by the system and drive the transmon from its ground state to its third excited state (|T3) in the main text). This 
photon-transmon interaction, not usually present in cQED devices, can, however, be achieved by using a nonlinear 
capacitative element. These capacitors have been realized and studied in multiple disciplines and with different 
implementations, such as using ferroelectric thin films or ceramics [4, 5], quantum wells in heterojunctions [6], MOS 
junctions [7, 8] and carbon nanotubes [9, 10]. 

We are interested in a nonlinear capacitor whose energy is a non-quadratic, symmetric and increasing function of 
the voltage for small fluctuations of V. In this way, the system will be confined into a region of small potentials [11]. 

Among the aforementioned nonlinear capacitors, the one based on carbon nanotubes seems to be the most suitable 
candidate. This capacitor [9] consists of two parallel plates with carbon nanotubes in between, placed perpendicular 
to the plates. The nonlinearity in the C(V) curve is due to the finiteness of the density of states (DOS) in the 
nanotubes [9, 10]. The capacitance of this device goes as C(V) ^ 00(1 -I- aU^), thus the energy spectrum of this 
capacitor as a function of the potential difference across the plates has the form 

E{V)^ dV'{V-V')C{V') =-CoiV^+aV'^), (S15) 

Jo 2 

where, to simplify the notation, the constant a has been redefined to absorb a numerical constant. 

Note that the carbon nanotubes do not constitute a key element in our proposal. Any other capacitative element, 
easier to fabricate and implement, with the same or similar energy spectrum will work as well. 


THE MODEL: DERIVATION 


The Lagrangian describing our circuit QED system is found in the same way as in circuit theory [12, 13]. The 
Lagrangian of the device presented in FIG. S4 is given hy L = Ls + Lt, where Ls is the Lagrangian of the SQUIDs, 


given by 


j- C2a , . TA \2 , G2sa , . . n2 , — ^2a\ 

Ls =—;^ {‘P2a -V2) H-((/3l - ^2a) + Ej2a COS - 

2 2 \ (/Jo y 

, ^26 , . TA \2 , (^286 , . . n2 , f — y^2b\ 

+ [^2b -V 2 ) H-— ((/Jl - ip2b) + Ej2bCOS - 

2 2 \ V’o / 

, C'sa / . TA ^2 , C'ssa / . . n2 , (‘fl — ^3a\ 

2 2 \ (/Jo / 

+ ^ {p3b - v^f + ((/Jl - (/Job)^ + Sjsb cos (——, (S16) 

2 2 \ (/Jo y 

where the dynamic variables (/Ji, (/J 2 q, (p 2 bj ‘P 3 a and ‘P 3 b are the fluxes defined at each node of the diagram in FIG. S4 
and (/Jo = h/{2e) is the flux quantum divided by 27r. For the transmon, the Lagrangian that describes its behavior is 
given by 

Lt =— {pt — Vif + — ((/j/ + ap\) + Ejt cos ■ (S17) 


Note that the expression in Eq. (S17) contains the energy spectrum of a carbon nanotube nonlinear capacitor [right- 
hand side in Eq. (S15)]. 

The Legendre transformation of the Lagrangian in Eq. (S16) and (S17), needed to study the quantum dynamics of 
the router, gives a function that is not easy to work with due to the presence of quadratic and cubic radicals. Instead, 
a Taylor expansion of the transformation is more suitable for our purposes. The expanded Hamiltonian is given by: 


-h 2p2aC2aV2 - C2saC2aVi 
2{C2a+C2sa) 

PSa + 2p3aC3aV3 " CssaCsgV^ 

2{C3a + C3sa) 
l3x‘^ P^x^ P^x^ 

1274 I 877 “ I 87 I 0 ~ 


Plb 


Plb 


‘2p2bC2bV2 — C2sbC2bV2 
2,{C2b + C2sb) 

^PsbCsbVs — CssbCsbV^ 


2(C'3b -I- Cssb) 

CiV,^ p fpi 

— - Ejt cos — 

2 VT’o 


- Ej 2 a cos 

- Ej3a COS 


/ Pi - P2a 

V 7>0 

/ Pi - P3a 

V PO 


Ej2b COS 


Ej3b COS 



(S18) 


For the sake of simplicity, the Hamiltonian in Eq. (S18) is expressed in terms of the fluxes, their conjugate momenta 
and the variable x, defined as 

X = pi + CiFi + —^f^(p2a + C2aV2) + „ {p2b + ^ 26 ^ 2 ) 

^2a ~r 02sa ^2b ~r 02sfc 

(7-3 (7-3 ^ 

+ ^ 7^ (P3a + G3a^3) + ^ + CafeHa), (S19) 

^ 3 a \ ^ 3 sa ^ 3 b “r ^ 3 sb 

with 


P = 6 Cta 


7 


Cl + Ct + 


C2aC2sa 
C2a + ^280 


C2bC2sb GaaGasa 
C2b + C2sb C3a + Gaso 


C3bC3sb 

Gab + Cssb 


(520) 

(521) 


In order to give rise to a density-density interaction between the different energy levels, the Hamiltonian is required 
to contain a term with PiP 2 a (see the next section in which the Hamiltonian is quantized). This expression appears 
after expanding the Hamiltonian in a Taylor series. In this same expansion a term containing p\Vi appears which 
describes the coupling between the third energy level of the transmon and the transmission lines. 

After substituting Eq. (S19) in Eq. (S18), many more interactions appear, but they can be neglected by imposing 
constraints on the parameters of the system. To start with, the coupling between the transmon and the outgoing 
transmission lines can be neglected by making G 28 a, G 28 b, Ga8a and Ga8b smaller than Gi, G 2 a, G 2 b, Gaa and Gab- In 
this way, the exchange interaction between the SQUID 2a and the third outgoing transmission line is also negligible. 
In order to see how the remaining couplings are enhanced or neglected we first have to quantize the Hamiltonian. 
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Csa ^ ^ Csb 

</?3a ifiSb 


Vi 


C2s 


Cssa i 'AEjSaX Cssb i >'EjSb'>' 


Cl 


<P1 


C2a i 


=F tEj2at :i(:Ej2bt 


‘P2a 


^2b 


C2b 


V 2 


Ejt 


H" 


Ct 


FIG. S4. Circuit QED device proposed to operate as a quantum router. This device is composed of four SQUIDs, each with 
a capacitor {C 2 sa, C 2 ab, Csaa, Csab) uud a pair of Josephson junctions (with energies Ej 2 a, etc.). It also contains a transmon 
qnbit, with a capacitor (Gt) and a pair of Josephson junctions (Ejt). These five elements are capacitively coupled to an 
incoming transmission line (Vi) and to two outgoing transmission lines (V 2 and V 3 ). The fluxes ifii, in the nodes of the circuit, 
are quantized variables (see the text around the figure, section The model: derivation, for the quantization of these variables) 
and describe the absorption and emission of incoming photons by the transmon and SQUIDs. In this proposal, two SQUIDs 
are placed in each branch for scalability purposes and robustness. Since the target photons can be in two states with different 
frequencies (for scaling purposes, for these will be the control photons in a subsequent router), two resonators (per branch) 
coupled to the same transmission line, with two distinct threshold energies for which they are detuned, are needed to absorb 
the target photons. 


Quantization of the Hamiltonian 


The Hamiltonian is quantized as usual by introducing the ladder operators as follows: 


Pi = 


n = 


ih 

2 lpq 

2I3Et 

3^^ 


G V'‘( c 


1/4 


(ai + a\). 


(522) 

(523) 


Here we have made use of the variable cjii = ipijipo, to simplify the derivations, where To = -^ and also Et = g~^- 

The operators oi and a\ are the usual ladder operators that annihilate and create an excitation in the transmon. For 
the SQUIDs, the quantized forms of the operators are 


P2a = - 


<t>2a — 


ih 

2 (fill 

I 


Ej2a 


C 2 


C2a + C2sa 


c- 


2sa 






1 / 2 ' 


1/2 


(a2a - 4a) 


Ej2a C2a + C2sa 




1/2 


(Oi + 4)- 


(524) 

(525) 
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Similar expressions can be derived for the other operators. With Eqns. (S22) - (S25) together with Ejt + Ej 2 a + 
Ej 2 b + Ej 3 a + Ejsb = 37^//3 = Ej and Ej 2 aE 2 a = {Ej - Ej 2 a + Ej 2 b + Ej 3 a + Ej 3 b) (for simplicity we have also 
introduced i? 2 a = SQUID-transmon exchairge interactions are strongly suppressed aird the expressions 

describing the energy levels of the system together with the density-density interactions between them are found [see 
Eq. (2) in the main text]. The SQUID-SQUID exchange interaction can be neglected since its interaction strength 
goes as (C' 2 sa/C' 2 a)^- Other similar interactions can be neglected by making the capacitances C 2 sa, C 28 b, Ossa and 
Cssb smaller than all the others and using the rotatirrg wave approximation. 

The interaction Hamiltonian [Eq. (3) in the main text], which couples the transmon and the SQUIDs with the 
incoming and outgoing transmission lines, is derived by imposing the following quantization of the potentials Ui, V 2 
and V 3 : 


= (S26) 

2(^0 

Here 6 ^, h\ are the ladder operators describing the oscillating modes in the i-th transmission line, and Ai are three 
constants that may be different. 

Let us look again the expressions in Eqns. (S18) and (S19). In the Hamiltonian [Eq. (S18)], there are only even 
powers of x. If we consider a single incoming photon, contained in Vi [see Eq. (S26)], given that the terms in the 
Hamiltonian containing a single power of Vi only contains odd powers of the momenta, the photon will excite only 
odd levels of the system. An even level can only be excited if an odd level is already occupied. This is shown in the 
Hamiltonian He in the main text [Eq. (3)]. 


Projection operators 


In order to simpliby both the notation and the calculations we have introduced three projection operators defined 


as 


t 2 t" 

“ti =g«l«l 


J.T2 


'3^/2 


aia[ = 


0 0 0 0 \ 
0000 
0001 
V 0 0 0 0 y 

/ 0 0 0 0 \ 
0001 
0000 
V 0 0 0 0 y 





/ 0 0 0 1 \ 
0000 
0000 
V 0 0 0 0 y 


(S27) 


(S28) 


(S29) 


to replace oi and its hermitian conjugate. The Hamiltonian in the main text is written in terms of these new operators 
(see the section Analysis). These operators have the property that acting on the ground state of the transmon 
creates an excitation in the i-th level, whereas the same operator acting on any other state where the transmon is 
already excited gives zero. Similarly, the operator aTi only annihilates an excitation of levels in the transmon, lowering 
it to its ground state. 

With these projection operators, the quantized Hamiltonian iLc, which describes the coupling between the trans¬ 
mission lines with the transmon and SQUIDs takes the simple form 


He = 



^ ) a^ 2 «Ti 6 i(p) + ( \ — -3 

V TrTTs/ VV tttti 


TTTTS 


a%3aT2hi{jp) 


{hiij)) + h2{p))) 


^3a 


/•^Ta 


+ 


^3b 




(blip) + b 3 {p))) + h.c. 


(S30) 
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The fourth term in the right-hand side, whose hermitian conjugate describes the process where the third excited level 
of the transmon decays into the second while emitting a photon, can be canceled by setting a as 


a = 




i2AhCt ■ 


(S31) 


In this way, the number of states the third level of the transmon can evolve into is reduced and, thus, the lifetime of 
the third level of the transmon is enhanced, increasing the transition matrix element from the ground state to |T3). 


Langevin equations of motion 

The Langevin equations of the ‘static’ variables of the system in terms of ^ 2 m and h'Hn are 


clti = — I iujTi + 
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(S32) 


(S33) 


(S34) 


where the index k in the sum runs over k £ {2a, 26,3a, 36} and three more equations for 26, 3a and 36, similar to 
Eq. (S34), have to be derived. Relaxation and dephasing ratios are then introduced into this set of Langevin equations 
by following Refs. [14, 16]. 


PHOTON MODES IN THE TRANSMISSION LINES 


In this section we briefly describe the meaning of the left/right-moving and even/odd modes, which are common in 
many circuit-QED works where the input-output formalism is used. For more details we refer the reader to Ref. [17], 
appendix A, where the following description is taken from. 

In case the photons in a transmission line are allowed to propagate in both directions, the photon modes are 
described by right-moving bosonic operators (r, rl) if they propagate to the right or by left-moving bosonic operators 
{I, /I) if they propagate to the left. In this case, the Hamiltonian that describes the transmission lines reads: 


Ht 



(S35) 


for each transmission line. The minus sign appears because the dispersion relation of the left-moving modes is the 
mirror image of the one of the right-moving modes. The remainder of the Hamiltonian given in the main text should, 
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Ct 

C2a 

C2b 

Caa 

Cab 

C-2sa 

C2sb 

Cssa 

Cash 

0.1 Cl 

1.2 Cl 

1.1 Cl 

Cl 

Cl 

0.1824 Cl 

0.1934 Cl 

0.1560 Cl 

0.1799 Cl 


TABLE S3. Proposal for a set of parameters that satisfy the requirements of a quantum router, given in terms of the capacitance 

_Q 2 

Cl = 10 ^F. With these values, the device can operate as described in the main text, forwarding photons into one of the 

outgoing transmission lines according to the state of the control photon. The rest of the parameters (Josephson energies, 
inductances, a) are found as functions of the capacitances during the derivation of the quantized Hamiltonian (see the second 
section of these supplementary notes) . 


in this case, contain the sum of operators r + l and +/f instead of b and for each transmission line, so the different 
components of the router interact equally with both the left and right-moving modes. 

From this new Hamiltonian, and following the same procedure as before (input/output formalism), similar equations 
of motion are found. Nevertheless, we can also do one more step and define two new operators: the even (6, b^) and 
odd (&, ¥) operators that describe what are called even and odd modes, defined as: 


, Tp -f l—p 

i _fp — l-p 

Ur) -" 


• 


(536) 

(537) 


These we then introduce in the Hamiltonian replacing the right and left-moving modes. One can check that the odd 
modes are not coupled to the system (they do not appear in the interacting part of the Hamiltonian), so they can be 
omitted in the Hamiltonian. The even modes are the ones used in the main text. 


NUMERICAL INTEGRATION OF THE SCATTERING AMPLITUDES 

In order to find the set of values for the capacitances —these are the only free parameters left— that best satisfy the 
requirements for an optimal quantum router it is convenient to maximize (or minimize) the scattering probabilities 
according to the operating requirements. Since an analytical solution has proved difficult to obtain, we have looked, 
by trial and error, for a set of numbers with which the router works as expected and found a set of capacitances 
that satisfy our requirements of energy spacing between the SQUIDs and transmon levels. Then, by maximizing the 
scattering amplitudes we have hne-tuned these values. Although a more systematic way of searching may result in a 
better solution, we have found that there exist at least one set of such parameters, shown in Table S3. The reason 
for this fine tuning is that, although the separation between the energy levels of the transmon is relatively large, the 
separation between the energy levels of the SQUIDs is very small (and the optimal operation of the quantum router 
sensitively depends on small variations of the values of the capacitances) 

With the values in Table S3, the energies are found to be around 10 GHz for the transmon energy levels {ujti = 
7.189 GHz and ujts = 21.92>GHz ) and around IGHz for the SQUIDs {uj 2 a = 1.028 GiJz, uj 2 b = 1.102GHz, 
W 3 o = 1.053 GiJz and uJzb = 1.186 GiJz). The coupling strength is of the order of 10~^ Hz. All these values are 
common in the related literature. Also, we used the ratio Ej/Et = 1296 which, although being larger than in other 
circuit-QED works, is sufficient to give rise to the nonlinear effects needed for the correct operation of the quantum 
router [3]. 

For the two-photon processes, we have made some assumptions: from the single-photon process, we have seen that 
the transmon absorbs the hrst incoming photon with a large probability, so for the two-photon processes we assumed 
that we sent the second photon while the transmon is already excited and that it does not decay emitting the control 
photon until the second photon has been fully transmitted. This implies that the lifetime of the transmon is much 
larger than that of the SQUIDs. This is a plausible assumption, given that the lifetimes of such cQED devices depend 
not only on the elements they are composed of but also on the geometry of the system (Purcell effect, see e.g. [18]). 

Finally, we comment on the possibility to use the quantum router in a larger network containing many quantum 
routers such that the transmitted photons of one router can be used as control photons in the next node of the 
network. This means that a photon transmitted by the router, which has enough energy to excite one of the SQUIDs, 
can also excite the transmon in the next router. Therefore, in a scalable network where all the routers are identical, 
the energy needed to excite the transmon (wti or wts) is the same as the energy needed to excite one of the SQUIDs 
(wa or iOb). In our (optimized) model, however, the energy of the hrst excited level of the transmon is one order of 
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magnitude larger than that of the SQUIDs with the current set of capacitances. A possible solution may be to add a 
device between the two consecutive routers that enlarges the energy of the first photon to arrive. Stimulated Raman 
Adiabatic Passage (STIRAP) processes may be a candidate for a device capable of transferring the population of one 
quantum state (with lower energy) to another (with higher energy). These processes have been studied in circuit 
QED systems in [19, 20]. 


REQUIREMENTS FOR AN EFFICIENT QUANTUM ROUTER FOR QUANTUM COMMUNICATIONS 

So far, different proposals for quantum switches or routers that may be implemented in a quantum network have 
been studied [21-25]. In order to assess the suitability of a given implementation for efficient use in a quantum 
communication network, we list here a set of requirements that we believe a quantum router must satisfy (part of 
these have been proposed elsewhere [24, 25]). The list below contains general requirements, which can be applied and 
checked for specific implementations. 

Imagine we have a network for quantum communication composed of quantum channels (which allow for the 
transmission of quantum information) and nodes (which contain quantum routers). Consider also that there is a 
signal to be sent and an element that controls the switches. 

1. The information in the signal to be transmitted and also the information encoded in the control elements must 
be stored in quantum objects (qubits). Otherwise we can not talk about a fully quantum-mechanical router [24]. 

2. The router has to be able to route the (quantum) signal into a coherent superposition of both output modes [24]. 

3. The signal information must remain undisturbed during the entire routing process [24]. The control information 
must also remain undisturbed, so that we can keep track of the signal and verify that it is not modified (due to 
e.g., entanglement with the control element). 

4. The router has to work without any need for postselection on the signal output [24] . 

5. To optimize the resources of the quantum network, only a single control qubit is required to direct each signal 
qubit [24]. Alternatively, the control information could also be stored in the signal qubit itself [25], but that 
may give the receiver information about the network which, in some applications, may not be a good idea. 

6 . With one single control qubit, as many as possible signal qubits should be forwarded. 

7. A router with n outputs must operate with the same efficiency (or speed, or with the same number of operations, 
etc.) for each output channel. 

8 . Scalability: the output of a router should be the input of the next one in a network. 

The quantum router that we propose satisfies all of these requirements except for the last one (although, as we 
discussed in the previous section, this problem can for example be amended by using STIRAP processes). We thus 
believe that it can become a good candidate for a quantum switch or router in a network for quantum communication. 


* salacadellans@gmail.com 
^ m.blaauboer@tudelft.nl 
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